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We analyze the three-mode correlation properties of the electromagnetic field in a optical para- 
metric oscillator below threshold. We employ a perturbative expansion of the Ito equations derived 
from the positive-P representation of the density matrix. Using the generalized Cauchy-Schwarz in- 
equality, we investigate the genuine quantum nature of the triple correlations between the interacting 
fields, since in this case continuous variable entanglement is not detected by the van Loock-Furusawa 
criterion [Phys. Rev. A 67, 0523f5 (2003)]. Although not being a necessary condition, these triple 
correlations are a sufficient evidence of tripartite entanglement. Of course, our characterization of 
the quantum correlations is applicable to non-Gaussian states, which we show to be the case of the 
optical parametric oscillator below threshold, provided nonlinear quantum fluctuations are properly 
taken into account. 



The nonlocal character of quantum mechanics has been 
widely exploited to develop quantum information pro- 
tocols [Ij] sometimes implemented in laboratories |2j], or 
even for commercialpurposes as has occurred with quan- 
tum cryptography '3*1 . Bipartite entangled states are the 
simplest resource for such protocols and can be generated 
by optical parametric oscillators (OPOs) \^±^. How- 
ever, more sophisticated protocols may require entangle- 
ment of many parts [SJ, |9|. The next difSculty is try- 
ing to find a source of tripartite entanglement. The first 
idea in this sense came from Greenberger, Home, and 
Zeilinger [lOj, suggesting the GHZ state, a discrete vari- 
able (spin variables) tripartite entangled state detectable 
by triple product spin measurements, but actually very 
difficult to implement. However the implementation of 
this idea by using continuous variables such as quadra- 
tures of optical fields, seems quite feasible with the avail- 
able homodyning techniques. Nowadays, there are a few 
devices that can generate multipartite entangled states, 
and the OPO has been explored for this task [11]. 

The nonlinear interaction between the OPO modes 
suggests the possibility of tripartite entanglement be- 
tween the three modes, and also the possibility of control- 
ling the bipartite correlations between signal and idler 
only by handling a third part, the pump mode. The 
quantum correlations between the different parts are sup- 
posed to change depending on statistics, spatial profile 
and other possible characteristics of the pump input, as 
observed in many experiments based on parametric down 
conversion 12|. Indeed, tripartite pump-signal-idler en- 



tanglement has been predicted in the intense regime, 
above the OPO threshold [13| and experimentally gener- 
ated [IJI . Tripartite correlations were analyzed according 
to the criteria proposed by van Loock and Furusawa|15j. 
based on pair correlations among the three modes. Also, 
there are other experiments that generate continuous- 
variable tripartite entangled states obtaine d by mixing 
squeezed vacua with linear optical elements |l6| . 

However, there has been little attention to tripartite 



correlations in the OPO operating below threshold, but 
it was in this context that triple quadrature correlations 
were proposed to be measured as means to compare clas- 
sical and quantum correlations 17|. The experimental 



difficulty resides in homodyning the weak down converted 
beams, which is not a problem for intense beams gener- 
ated above threshold, since a smart technique of self- 
homodyning 1J| can be used. It is also important to 



notice that it was in the below threshold regime where 
squeezed vacuum was first produced, and EPR correla- 
tions were first detected [ISl ll9| . 

In this work, we calculate, using a perturbative ap- 
proach, all triple quadrature correlations between the 
interacting fields in the OPO below threshold. These 
triple correlations should vanish for Gaussian quantum 
states, whether they are entangled or not, so that the 
van Loock and Furusawa|l5| entanglement criterion is 
well suited. This is the case, for example, of the OPO 
above threshold luj. However, we shall demonstrate 



that triple quadrature correlations are predicted for be- 
low threshold operation when the perturbative approach 
is taken beyond the linear approximation. While evi- 
dencing non Gaussian behavior, these triple correlations 
are also sufficient conditions for tripartite entanglement. 
We also show that, in this case, the van Loock-Furusawa 
criterion fails to evidence entanglement. Due to this in- 
conclusive condition, and the lack of a criterion for this 
case, we also show that certain correlations violate the 
tripartite Cauchy-Schwarz inequality expected for classi- 
cal c-number variables. 

The OPO can be modeled by three modes coupled 
through a nonlinear crystal inside a lossy Fabry-Perot 
cavity, described by the follo wing master equation in the 
Born-Markov approximation [20] : 
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The parameter E represents the external pump field at 
frequency loq. The pump, signal and idler fields are rep- 
resented by the annihilation operators oqi ^1 a-nd 02, 
respectively, and satisfy the frequency match condition 
LjQ — uji + UJ2. X is the nonlinear coupling and 7^ 
(i = 0,1,2) are amplitude damping rates. Now, the 
density matrix equation of motion can be transformed 
into c-number Fokker-Planck or stochastic equations by 
means of operator representation theory. 

The quantum evolution of the system can be described 
by the generalized P-representation [21|, by expanding 
the density matrix in a off diagonal coherent state basis, 
defined as 
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where a = (ao,ai,a2) and a^ = (qq ,Q!i ,0:2 ) are com- 
plex variables that run independently over all complex 
plane. The function P{a,a~^) can be understood as a 
positive distribution in a double phase space that satisfies 
a Fokker-Planck equation. The Ito stochastic equations 
derived from the Fokker-Planck equation, after trans- 
forming to the rotating-wave frame, are 



dao = 

duQ — 

duj — 

dat = 



{E - 7oao - XCH"2) dt 

(E* - joa^ - X'^iO^t) dt 

{-IjOtj + Xftfe ao) dt + {xaaf^ dWj 

(-7,a+ + x^kat) dt + (xa+) ^'^ dW+ , (4) 



where j = 1,2; k = 1,2 and j ^ k. The Wiener incre- 
ments satisfy {dWidW2) = {dW^ dW}) = dt , and ah 
other correlations vanish. 

In order to analyze the correlations between the modes, 
we proceed by applying the perturbation theory for the 
three modes below threshold, where this approach is 
valid. Around the threshold, the quantum fluctuations 
become huge and break the perturbative analysis. Fol- 
lowing the physical situation found in most experiments, 
we shall assume a common damping rate for the down 
converted fields: 71 = 72 = 7. It is also useful to define a 
dimensionless coupling constant g = x/ilV^lr), where 
7r = 70/7- In order to employ a perturbative approach 
and to decouple the dynamical equations, it is convenient 
to define the following normalized variables [2 2] : 
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Note that in positive-P representation the quadrature 
variables x and y are not real and their definitions are ac- 
companied by the corresponding phase space conjugates 
x'^ and y^ . We now expand these variables in terms of 
the dimensionless coupling constant |23l.l24|: 
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and analogous expansions for y^, x, y, a;+ and y+. 

Below threshold, the zero-order steady state solution 
for the quadratures are given by 
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where fi — x-E^/(77o) is the dimensionless pump parame- 
ter. These are the mean values of the fields below thresh- 
old. As expected, at this regime all mean field quadra- 
tures are zero, since only fluorescence is present. The 
pumped quadrature is the only one having a macroscopic 
value. The first quantum correction comes from the next 
order contribution: 
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where, {Ut)i,+ {t')) = {iy{t)S,y+{t')) = 6{t - t') . Note 
that we kept the second order perturbation term for 
the pumped mode variables to account for pump deple- 
tion, since their first order corrections vanish. Therefore, 
pump depletion is a second order effect below threshold. 
It is also interesting to note that the second order cor- 
rections in the quadrature variables of the pumped mode 
depend on the product of two correlated Gaussian pro- 
cesses. This implies that double correlations between the 
pump and any one of the down converted modes van- 
ish. As a result, a criterion based on pair correlations for 
tripartite entanglement, like the one by van Loock and 
Furusawa [l5|, fails to detect the existing triple correla- 
tions which are quantum mechanical in essence, as we 
shall see shortly. Moreover, the second order correction 



to the pump fluctuations are not Gaussian, since it can 
be checked that its odd moments do not vanish. 

With the presented solutions, we can calculate any 
kind of correlations between the mode variables. Cor- 
relations between two quadratures generate the expected 
two- mode squeezed state for yy^ . Actually, the four vari- 
ables of the down converted modes can be combined to 
form EPR variables. An important observation is that 
there are no two-mode correlations between any variable 
of the pumped mode and the down converted ones, since 
the pumped mode depend on the product of two Gaus- 
sian noises. Due to this property, the three-mode corre- 
lations unveil most of the important result addressed by 
this paper. 

The results for the steady state time ordered triple cor- 
relations calculated in normal order, is presented below. 
Actually these results are also the output cavity fields, 
since in the positive-P representation the incoming field 
(vacuum) is not correlated with the intracavity modes. 
The non vanishing triple correlations up to fourth order 
in the coupling constant are: 



(AxAx+Axo) 

{AyAy+Axo) 

(AyAx+Ayo) 
(AxAy+Ayo) 



fJ- 



1-M 

-fr 



1 + M 



9 



2 r 
4 / M 






7^ + 2 (1 + Ai) 

,2 



= 9 



^^ 



2 



1-M 



7r 



A* 



i-ny \2 + j, 



Ir 



(9) 



where Au = m — {u). One can easily show that for any 
kind of separable state of the form: "Y^PiPk ® Pirm where 
{k,l,m} is any permutation of {0, 1,2}, at least one of 
these triple correlations should vanish identically. There- 
fore, the set of triple correlations given by eq.® is a suf- 
ficient condition for tripartite entanglement. However, it 
does not exclude classical correlations [25|, so that fur- 
ther analysis is required to evidence genuine quantum 
behavior. 

Let us now turn to a reliable classicality criterion based 
on Cauchy-Schwarz inequalities. First, we calculate the 
average product of the down converted squared fiuctua- 
tions: 

{{Ax^ + Ay')[{Ax+r + {Ay+n^ 
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The fluctuations on xq are given by 
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while {(Ayo) ) is negligible. Now we establish a relation 
between correlations © and intensity correlations. To 
begin with, we assume that for any three complex num- 
bers the following relation holds 



\a1a2 + A*aoP) > , 
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so that 



(a*aia2a2)-t-|Ap(aoao)+A(aiQ;2ao)+A*(a*a2ao) > . 

(13) 
This expression is minimized for A = {a*a2aQ)/{aQao), 
giving rise to the following inequality: 



(a*aia2a2)(aoao) > |(aia2ao)|^ 



(14) 



In the same way, similar inequalities can also be derived 
by permutation of the indices in the complex variables 
aj. After a straightforward algebra, it is easy to show 
that the triple correlation appearing in the right hand 
side of inequality P^ is simply the squared sum of all 
triple correlations given by eq.®, taking the flrst one 
((Ax Aai+Axo)) with an inverted sign. Also, the left 
hand side of the inequality is simply the product of the 
averages given by egs. fTU)) and (|lip . 

Using the results given by eqs.®, (|10l) . and (ITT|) on 
both sides of the Cauchy-Schwarz inequality ^^, we can 
verify its violation for the three-mode state produced by 
the OPO operating below threshold, which means that 
these triple correlations are essentially quantum mechan- 
ical. The reliability of the analytical calculations can be 
checked by numerical simulation of the full set of stochas- 
tic equations ([1]), without any perturbative approx- 
imation. Fig.(l) compares both sides of inequality (fT4| . 
showing a clear violation of the inequality for large 7^.. 
The numerical results are in very good agreement with 
the analytical calculations. In our numerical calculations 
we used g = 0.0071 {x = 0.01) whose order of magnitude 
corresponds to most nonlinear media used for paramet- 
ric down conversion. Note that the violation of Cauchy- 
Schwarz inequality occurs with very small correlation val- 
ues (~ 10""'). This is due to the rather small value of the 
coupling constant. Since these correlations scale as g^, 
this result can be substantially improved with new mate- 
rials providing sgnificant values of the nonlinear suscep- 
tibility. Indeed, we checked this numerically. 

In conclusion, we investigated the triple correlations 
produced in an OPO below threshold, evidencing tripar- 
tite entanglement not detected using the well known van 
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Figure 1: Violation of Cauchy-Schwarz inequality for an OPO 
operating below threshold. The figures show simulations us- 
ing ^ = 0.7 and g — 0.0071. (a) Bottom curve (black) rep- 
resents the right hand side of inequality ([14]), and top curve 
(red online) represents the left hand side for jr ~ 0.01. The 
Cauchy-Schwarz inequality is not violated, (b) Top curve 
(black) represents the right hand side of inequality (|14p . and 
bottom curve (red online) represents the left hand side for 
7r = 100. Violation of the Cauchy-Schwarz inequality is clear. 



Loock-Furusawa criterion. In order to circumvent this 
limitation, we employed a non classicality criterion based 
on general Cauchy-Scliwartz inequalities for third order 
moments. This approach was applicable due to the non 
Gaussian character of the quantum fluctuations. More- 
over, these triple correlations rule out separable states, 
being therefore a sufficient condition for entanglement. 
Although these inequalities have already been used for 
discrete variables [10], this is the first time, to our knowl- 
edge, they are used to characterize non classical triple 
correlations in continuous variables. Whether these triple 
correlations can be useful for quantum information pro- 
tocols remains to be investigated. 
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